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Abstract 

This paper focuses on the study of certain classes of Boolean functions that have ap- 
peared in several different contexts. Nested canalyzing functions have been studied 
recently in the context of Boolean network models of gene regulatory networks. In 
the same context, polynomial functions over finite fields have been used to develop 
network inference methods for gene regulatory networks. Finally, unate cascade 
functions have been studied in the design of logic circuits and binary decision di- 
agrams. This paper shows that the class of nested canalyzing functions is equal to 
that of unate cascade functions. Furthermore, it provides a description of nested 
canalyzing functions as a certain type of Boolean polynomial function. Using the 
polynomial framework one can show that the class of nested canalyzing functions, 
or, equivalently, the class of unate cascade functions, forms an algebraic variety 
which makes their analysis amenable to the use of techniques from algebraic geome- 
try and computational algebra. As a corollary of the functional equivalence derived 
here, a formula in the literature for the number of unate cascade functions provides 
such a formula for the number of nested canalyzing functions. 
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1 Introduction 



Canalyzing functions were introduced by Kauffman 1(| as appropriate rules in Boolean net- 
work models of gene regulatory networks. The definition is reminiscent of the concept of "canal- 
isation" introduced by the geneticist Waddington |22J to represent the ability of a genotype to 
produce the same phenotype regardless of environmental variability. Canalyzing functions are 
known to have other important applications in physics, engineering and biology. They have 
been used to study the convergence behavior of a class of nonlinear digital filters, called stack 
filters, which have applications in image and video processing (5); [23; 24 1. Canaly zing functions 



also play an important role in the study of random Boolean networks |10|; [13|; [I9j; [20J , have 



been used extensively as models for dynamical systems as varied as gene regulatory networks 
lol |. evolution, 20] and chaos [13]. One important characteristic of canalyzing functions is 
that they exhibit a stabilizing effect on the dynamics of a system. For example, in [5J, it is 
shown that stack filters which are defined by canalyzing functions converge to a fixed point 



called a root signal after a finite number of passes. Moreira and Amaral [15|], showed that the 



dynamics of a Boolean network which operates according to canalyzing rules is robust with 
regard to small perturbations. 



A special type of canalyzing function, so-called nested canalyzing functions (NCFs) were intro- 
duced recently in Q , and it was shown in j^] that Boolean networks made from such functions 
show stable dynamic behavior and might be a good class of functions to express regulatory 
relationships in biochemical networks. Little is known about this class of functions, however. 
For instance, there is no known formula for the number of nested canalyzing functions in a 
given number of variables. 



Another field in which special families of Boolean functions have been studied extensively is 
the theory of computing, in particular the design of efficient logical switching circuits. Since 
the 1970s, several families of Boolean functions have been investigated for use in circuit design. 
For instance, the family of fanout-free functions has been studied extensively, as well as the 
family of cascade functions. A subclass of these are the unate cascade functions see, e.g., 



141 ; |16| |. which we focus on here. It turns out that this class of functions has some very useful 
properties. For instance, it was shown recently that the class of unate cascade functions is 
precisely the class of Boolean functions that have good properties as binary decision diagrams. 
In particular, the unate cascade functions (on n variables) are precisely those functions whose 

binary decision diagrams have the smallest average path length (2 — 7^— [) among all Boolean 
functions of n variables. 



The notion of average path length is one cost measure for binary decision trees, which measures 
the average number of steps to evaluate the function on which the tree is based. One way of 
assessing the relative efficacy of classes of Boolean function for logic circuit or binary decision 
tree design is to look at the number of different circuits or trees that can be realized with a 
particular class. That is, one would like to count the number of functions in a given class. This 
has led to a formula for the number of unate cascade functions One of the results in this 
paper shows that the classes of unate cascade functions and nested canalyzing functions are 
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identical (as classes of functions rather than as classes of logical expressions). As a result of the 
equivalence we will establish, this formula then also counts the number of nested canalyzing 
functions. 



A third framework for studying Boolean functions, in the context of models for biochemical 



networks, was introduced in [111 ] . There, a new method to reverse engineer gene regulatory 
networks from experimental data was proposed. The proposed modeling framework is that 
of time-discrete deterministic dynamical systems with a finite set of states for each of the 
variables. The number of states is chosen so as to support the structure of a finite field. One 
consequence is that each of the state transition functions can be represented by a polynomial 
function with coefficients in the finite field, thereby making available powerful computational 
tools from polynomial algebra. This class of dynamical systems in particular includes Boolean 
networks, when network nodes take on two states. It is straightforward to translate Boolean 
functions into polynomial form, with multiplication corresponding to AND, addition to XOR, 
and addition of the constant 1 to negation. In this paper we provide a characterization of those 
polynomial functions over the field with two elements that correspond to nested canalyzing 
(and, therefore, unate cascade) functions. Using a parameterized polynomial representation, 
one can characterize the parameter set in terms of a well-understood mathematical object, a 
common method in mathematics. This is done using the concepts and language from algebraic 
geometry. To be precise, we describe the parameter set as an algebraic variety, that is a set 
of points in an affine space that represents the set of solutions of a system of polynomial 
equations. This algebraic variety turns out to have special structure that can be used to study 
the class of nested canalyzing functions as a rich mathematical object. 



2 Boolean Nested Canalyzing and unate cascade Functions are equivalent 



2. 1 Boolean Nested Canalyzing Functions 



Boolean nested canalyzing functions were introduced recently in [8|, and it was shown in 
[qI that Boolean networks made from such functions show stable dynamic behavior. In this 
section we show that the set of Boolean nested canalyzing functions is equivalent to the set 
of unate cascade functions that has been studied before in the engineering and computer 
science literature. In particular, this equivalence provides a formula for the number of nested 
canalyzing functions in a given number of variables. 

We begin by defining the canalyzing property. 

Definition 2.1 A Boolean function f(x±, . . . , x n ) is canalyzing if there exists an index i and 
a Boolean value a for Xi such that f(x\, . . . , a, Xj+i, . . . ,x n ) = b is constant. That is, 
the variable x^, when given the canalyzing value a, determines the value of the function f , 
regardless of the other inputs. The output value b is called the canalyzed value. 
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Throughout this paper, we use the Boolean functions AND(x, y) = x Ay, OR(x,y) = x V y 
and NOT(a;) = x. 



Example 2.2 The junction AND(x, y) = x Ay is a canalyzing function in the variable x with 
canalyzing value and canalyzed value 0. The function XOR(x, y) := (x V y) A (x A y) is not 
canalyzing in either variable. 

Nested canalyzing functions are a natural specialization of canalyzing functions. They arise 
from the question of what happens when the function does not get the canalyzing value as input 
but instead has to rely on its other inputs. Throughout this paper, when we refer to a function 
of n variables, we mean that / depends on all n variables. That is, for 1 < i < n, there exists 
(a 1 , ...,a n )eF% such that f(a 1 , . . . , a^, a u a i+1 , ...,a n )^ f{a 1 , . . . , Oj_i, a7, a i+1 , . . .,a n ). 

Definition 2.3 Let f be a Boolean function in n variables. 

• Let a be a permutation on {1, . . . ,n}. The function f is nested canalyzing function (NCF) 
in the variable order x a ^, . . . , x„( n ) with canalyzing input values and canalyzed 

output values b\, . . . , b n , respectively, if it can be represented in the form 

h if = ai, 

b 2 if x a (i) ^ ai and x a ( 2 ) = a 2, 
b 3 if x ai i) ^ ai and x a ( 2 ) ^ a 2 and x a ( 3 ) = a 3 , 



f{x x ,x 2 , ...,x n ) = < 



(2.1) 



b n if x a (i) ^ ai and • • • and x a (n-i) ^ On-i and x a ( n ) = a n , 
b n if x a (i) ^ ai and • • • and x^ n ) ^ a n . 

• The function f is nested canalyzing if f is nested canalyzing in the variable order x^x), ■ ■ ■ , x<r{n) 
for some permutation a. 

Example 2.4 The function f(x,y,z) = x Ay A z is nested canalyzing in the variable order 
x,y,z with canalyzing values 0,1,0 and canalyzed values 0,0,0, respectively. However, the func- 
tion f(x,y,z,w) = x Ay A XOR(z,mj) is not nested canalyzing because if x = 1 and y = 1, 
then the value of the function is not constant for any input values for either z or w. 

The following lemma follows directly from the definition above. 

Lemma 2.5 A Boolean function f on n variables is nested canalyzing in the variable order 
x a m , . . . , x c ( n \ with canalyzing input values and canalyzed output values b n , 

respectively if and only if f depends on all n variables and, for all 1 < i < n, 

/(ci, . . . ,c„) = bi, 

where (ci, . . . , c„) G F2 such that c a ^ = and, for 1 < j < i, c a (j) = a]. 

The next lemma gives the functional form of a nested canalyzing function. To simplify nota- 
tion, we will use the following notational convention. Let a = 0, 1. Then x + a will denote x 
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if a = and x if a = 1. 
Lemma 2.6 Let 

g(xi, . . . , x n ) = (z CT (i)+ai+&i)Oi((z CT ( 2 )+a 2 +& 2 )<>2(- ■ ■ ((a; ff ( n _i)+a n _i+6 )l -i)^ ri -i(i ff ( ri )+a n +& n )) ■ • •), 

(2.2) 

u>/iere 

[v, ifbi = l; 

Oi=\ (2.3) 
[A, ^ = 0, 

and aj,&j G {0,1} /or all 1 < i < n. Then g is nested canalyzing in the variable order 
. . . , x a (n) with canalyzing input values a%, . . . , a n and canalyzed output values b\, . . . , b n , 
respectively. Furthermore, any nested canalyzing function can be represented in the form (12.21) . 

PROOF. It is clear that g depends on all variables x 1; . . . , x n . Let g n = £ CT („) + a n + b n and, 
for 1 < i < n, let 

9i = +di + bi)Oig i+1 . 
Then g = g x = (x ct( i) + a x + 6i)Oi^2- 

If x a (i) = ai, then (x CT (i) + a x + 6i)Oio 2 = &iOig 2 — W, by equation (12. 3p . For 1 < i < n — 1, 
suppose x a (j) = Hj for j < i and = a^. Now, for all j < i, we have bjOjgj+i = o J+1 and 
bi^igi+i = h. Thus, by equation (I2.3p . we get 

6iOi(6 2 2 (. . . (biOig i+1 ) ■■■)) = bi<>ig l+1 = b { . 

Hence g is nested canalyzing, with the as canalyzing values and the 6j as canalyzed values. 
It is left to show that any nested canalyzing function can be represented in the form (12. 2p . Let 
/ be a nested canalyzing function in the variable order a^i), . . . , £ CT (n) with canalyzing input 
values ai, . . . , a n and canalyzed output values b\, . . . , b n , respectively. By Lemma 12.51 and the 
above, it is clear that /(ci, . . . , c n ) = p(ci, . . . , c n ) for all (ci, . . . , c n ) G F^. Thus / = g as 
functions and hence / can be represented in the form (I2.2p . 



2.2 NCFs are Unate Cascade Functions 



We next show that Boolean NCFs are equivalent to unate cascade functions. Unate cascade 
functions have been defined and studied [l4|; llo] as a special class of fanout-free functions 
which are used in the design and synthesis of logic circuits and switching theory 

Definition 2.7 A Boolean function f is a unate cascade function if it can be represented as 

f(x!, x 2 ,..., x n ) = <(i ) 1 « (2) 2 (. . . (a£ (n _i)O n _ia£ (n) )) . . .), (2.4) 

where a is a permutation on {1, ... , n}, x* is either x or x + 1 and O, is either the OR (V ) 
or AND (A) Boolean operator. 
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Theorem 2.8 A Boolean function is nested canalyzing if and only if it is a unate cascade 
function. 



PROOF. Let / be a unate cascade function in the form (12.41) . Let a n and b n be such that 
x *u(n) = x o-(n) + a n + b n and, for 1 < % < n, let 

fl, ifO< = V; 
bi={ ' (2.5) 
[o, ifO i = A, 

and let a* G {0, 1} such that x*^\ = x a ^) + + b{. That is, 

a i = X l(i) + X <?(i) + b i- (2.6) 

Then 

f(x u ...,x n ) = {x a{l) +ai+b 1 )0 1 ((x a{ 2)+a 2 +b 2 )0 2 {- ■ ■ ((x CT ( n _i ) +a ri _i+6 n _i)O n _i(a; (T ( ri )+a n +6 n )) ■ ■ ■ ), 
which is nested canalyzing by Lemma 12.61 

Conversely, let / be a nested canalyzing function of the form (12. ip . By Lemma [2761 and equation 
(12. 2p . / can be represented in the form (12.41) where x*/* = x a ^ + + hi for all 1 < % < n. 
Thus / is unate cascade. 

Remark 2.9 The second sentence in the proof above implies the nested canalyzing function 
with canalyzing input (ai, . . . ,a n ) and canalyzed output (bi,...,b n ) is also nested canalyz- 
ing in the same variable order with canalyzing input (ai, . . . ,a n _i,a^) and canalyzing output 

6„_i,6„)- 

Remark 2.10 The theorem above provides natural equivalence between the class of nested 
canalyzing functions and that of unate cascade functions. Namely, for a given unate cascade 
function in the form ft2.4\ ), using \2. 6]) and \2. 5\) we can explicitly define the canalyzing input 
values ai and canalyzed output values bi. On the other hand, any nested canalyzing function 
in the form (OOP can be presented as a unate cascade function in the form \2.1$ where x*^ = 
x a(i) + a i + bi and Oj as in A2.3\) for all 1 < i < n. 

Using the theorem and remark above, it is now possible to translate results about one type 
of function into results about the other type. We point out one such example. Several papers 
have been dedicated to counting the number of certain fanout-free functions including the 
unate cascade functions [H; 0; @; [l7; 18]. On the other hand, Just et. al. 0, gave a formula for 
the number of canalyzing functions. 



Bender and Butler [2| and Sasao and Kinoshita [18( independently found the number of unate 
cascade functions, among other fanout-free functions. As an immediate corollary of Theorem 
12.81 we therefore know the number of NCFs in n variables, for a given value of n. We use the 



recursive formula found by Sasao and Kinoshita [18|, in the following corollary. 
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Corollary 2.11 The number of NCFs in n variables, denoted by NCF(n), is given by 

NCF(n) =2-E(n), 

where 



E(l) = l, £(2) =4, 



and, for n > 3, 



n-l 



E{n) = £ 



r=2 



n 

>r — 1. 



2 1 '- 1 .E(n-r + l) + T. 



For example, the number of Boolean NCFs (unate cascade functions) on n variables for n < 8 
is given by Table [Tlwhich is part of the tables given by Sasao and Kinoshita 18] and also 
Bender and Butler [2J. 



Table 1 

The number of NCFs on n < 8 variables 



n 


1 


2 


3 


4 


5 


6 


7 


8 


NCF(n) 


2 


8 


64 


736 


10,624 


183,936 


3,715,072 


85,755,392 



Some interesting facts derive from the equivalence of NCFs with unate cascade functions. 
Sasao and Kinoshita 0, Lemma 4.1], found that unate cascade functions are equivalent to 
fanout-free threshold functions. Thus NCFs are a special class of threshold functions. It is also 
interesting to note that among switching networks, the unate cascade functions are precisely 
those with the smallest average path length, as shown by Butler et. al. [3], which makes them 
efficient in logic circuit design. 



3 Nested canalyzing functions as polynomial functions 



Wanting to compute the total number of Boolean functions of a particular type, e.g. nested 
canalyzing or unate cascade functions, is one example of the need to study the totality of 
such functions. Few tools other than elementary combinatorics are available for this purpose, 
however. In this section, we propose an alternative approach to Boolean functions which 
provides a whole new set of mathematical tools and results. We will view Boolean functions 
{0, l} n — ► {0, 1} as polynomial functions / : FJ? — > F2, where F2 denotes the field with 
two elements. It is well-known that any function / : — ► F2 can be represented as a 
polynomial function 12j]. If we require that every variable appear with exponent 1, then 
this representation is unique. For Boolean functions, this representation is straightforward to 
construct by observing that AND(x,y) = x Ay = xy, OR(x,y) — x\/y — x + y + xy, and 
NOT(x) = x = x + 1. Conversely, replacing multiplication by AND and addition by the XOR 
function, we can translate any polynomial function into a Boolean function. In particular, 
this shows that any binary function on n variables can be represented as a Boolean function. 
While this seems like a simple change of language, it has the profound effect of placing the 
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study of Boolean functions into the fields of algebra and algebraic geometry, which have a rich 
body of results and algorithms available. 

The goal of this section is to formalize this equivalence and to characterize those polynomial 
functions that represent nested canalyzing functions. The characterization will be expressed 
as a parametrization, with the set of parameters taken from an algebraic variety. Algebraic 
geometry has many tools to study varieties, an approach that will be pursued elsewhere. 

3. 1 Polynomial form of nested canalyzing functions 

We derive a polynomial representation of the class of Boolean nested canalyzing functions 
which we then use to identify necessary and sufficient relations among their coefficients. 

Any Boolean function in n variables is a map / : {0, 1}" — ► {0, 1}. The set of all such maps, 
denoted by B n , can be given the algebraic structure of a ring with the boolean operators XOR 
for addition and the conjunction AND for multiplication. 

Consider the polynomial ring F2[xi, . . . , x n ] over the field F2 := {0, 1} with two elements. Let 
/ be the ideal generated by the polynomials x\ — Xi for alH = 1, . . . ,n. (That is, / consists 
of all linear combinations of these polynomials with arbitrary polynomials as coefficients.) 
For any Boolean function / G B n , there is a unique polynomial g G F 2 [x±, . . . ,x n ] such that 
g(ai, . . . , a n ) = f(a±, . . . , a n ) for all (a±, . . . , a n ) G F2 and such that the degree of each variable 
appearing in g is equal to 1. Namely 

n 

g(x u ...,x n ) = J2 f(a 1 ,...,a n )Y[(l-(x i -a i )). (3.1) 

(oi,...,a„)6F5 i=l 

And it is straightforward to show that this equivalence extends to a ring isomorphism 

R := ¥ 2 [xi, . . . ,x n }/I = B n . 
From now on we will not distinguish between these two rings. 

Next we present and study the set of all Boolean nested canalyzing functions as a subset of 
the ring R of all Boolean polynomial functions. 

The following theorem gives the polynomial form for canalyzing and nested canalyzing func- 
tions. 

Theorem 3.1 Let f be a function in R. Then 

(1) The function f is canalyzing in the variable X{, for some i, 1 < i < n, with canalyzing 
input value ai and canalyzed output value hi, if and only if 

f(xi, x 2 ,...,Xi,...,x n ) = (xi - a i )g(x 1 ,x 2 , . . . , x u . . . , x n ) + (3.2) 
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(2) The function f is nested canalyzing in the order xi,X2, ■ ■ ■ ,x n , with canalyzing values 
cii and corresponding canalyzed values b^, 1 < i < n, if and only if it has the polynomial 
form 



f(xi,x 2 , ■ . . , x n ) = (xi - ai)[(a; 2 - a 2 )[. . . [(x„_i - a n _i)[(x n - a n ) 

+ipn - 6 n -i)] + (K-i ~ 6«- 2 )] • • •] + (62 - 61)] + h 



(3.3) 



or, equivalently, 



n-l 



f(xi, x 2 , ■ ■ ■ , x n ) = Y[(Xi - Hi) + 

i=l 3=1 



n-j 

(b n - j+ i - b n -j) J[ (xi - a,i) 
i=i 



+ 61. (3.4) 



PROOF. 



(1) It is easy to see that if Xi — then the output is 6j, no matter what the values of the 
other variables are. Conversely, if / is canalyzing with input a« and output bi, then f — bi, 
as a polynomial in Xi has aj as a root, hence is divisible by X{ — Oj. This proves the first 
claim. 

(2) Let / be a nested canalyzing function as in Definition 12.31 and let 

g(xt, ...,x n ) = (xx- ai)[(x 2 - o 2 )[. • • [{x n -i ~ a n -i)[{x n - a n ) 

+ ipn - &n-l)] + {bn-1 ~ b n - 2 )\ ■ ■ •] + (&2 ~ &l)] + &1- 

We will show that ^ is the unique polynomial presentation of /, as in equation 13. II Since 
the degree of each variable in g is equal to one, we only need to show that g(ci, . . . , c n ) = 
/(ci, . . . , Cn) for all (ci, . . . , c n ) G F 2 . 

Clearly, if c\ = a±, then gf(ci, . . . , c„) = 61 . If c\ ^ a\ and c 2 = a 2 , then (ci— ai) = 1 and 
g(ci, . . . , c n ) = b 2 . If Ci ^ ai, c 2 ^ a 2 and c 3 = a 3 , then p(ci, . . . , c n ) = 63. We continue 
until we have c, 7^ for all 1 < i < n and c n = a n , in which case we get #(ci, . . . , Cn) — b n . 
If Ci 7^ ai for all i, then (q — Oj) = 1 for all i and hence g(ci, . . . ,c n ) = 1 + b n . Thus is 
the unique polynomial representation of /. 



3.2 A Parametrization of NCFs 



Our next goal is to derive a criterion as to when a given Boolean polynomial in n variables 
is a nested canalyzing function. The criterion will be given in terms of a parametrization 
of such polynomials corresponding to points in the affine space F^" that satisfy a certain 
collection of polynomial equations. Such a set is by definition an algebraic variety, in the 
language of algebraic geometry. This parametrization describes the entire space of nested 
canalyzing functions as a geometric object, whose properties can then be studied with the 
tools of algebraic geometry. 
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Recall that the ring of Boolean functions is isomorphic to the quotient ring R = ¥2[xi, . . . , x n ]/I, 
where / = {xj—Xi : 1 < i < n). Therefore, the terms of a Boolean polynomial consist of square- 
free monomials. Thus, we can uniquely index monomials by the subsets of [n] := {1, . . . ,n} 
corresponding to the variables appearing in the monomial, so that we can write the elements 
of R as 

R = {H c sU x i ■ cseF 2 }. (3.5) 

SC[n] ieS 

As a vector space over F 2 , R is isomorphic to Flf* via the correspondence 

R3 X^n^ — ► (c0, • • • , qn]) e Iff, (3.6) 

SC[n] ies 

for a given fixed total ordering of all square-free monomials. That is, a polynomial function 
corresponds to the vector of coefficients of the monomial summands. In this section we identify 
the set of nested canalyzing functions in R with a subset V ncf of F| n by imposing relations on 
the coordinates of its elements. 

Let S be any subset of [n]. We introduce a new term called the completion of S. 

Definition 3.2 Let S be a a non-empty set whose highest element is r$- The completion of 
S, which we denote by [rs], is the set [rs] ■= {1,2, . . . , rs}. For S = 0, let [r$] := 0. 

The main result of this section is the following theorem. 

Theorem 3.3 Let f be a Boolean polynomial in n variables, given by 

f(x h x 2 ,. . . ,x n ) = °s H x i (3-7) 

SC[n] ieS 

The polynomial f is a nested canalyzing function in the order 21,22, ■ ■ ■ ,x n if and only if 
C[ n i = 1, and for any subset SC. [n], 

°s = c [r s ] II C N\R- ( 3 - 8 ) 

ie[r s ]\s 

PROOF. First assume that the polynomial / is a Boolean nested canalyzing function in 
the order 21,22, . . . ,x n , with canalyzing input values a* and corresponding canalyzed output 
values hi, 1 < % < n. Then, by part 2 of Theorem 13.11 / has the form ( 13.4ft which can be 
expanded as 

n-X 

f(xi, . . . ,x n ) = Tl x i II a i + J2(b n -j+i-b n -j)( Tl x i II a i) + h- (3.9) 

SC[n]ieS le[n]\S 3=1 SC{n-j]iGS i e [„_j]\5 

We now equate corresponding coefficients in equations ( 13.71) and ( \3.9\) . First let S — [n]. Then, 
clearly, q rs ] = 1. Next, consider subscripts of the form S = [n]\{i},i 7^ n, that is, coefficients 
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of monomials of total degree n — 1 which contain x n . It is clear from equation (13.91) that x n 
only appears in the first summand and hence, for 1 < i < n — 1, 

c [n]\{i} — O^i — c [r s ]C[n]\{i}- (3.10) 

It is easy to check that equation (13.81) holds for any set 5* C [n] such that S = [rs] ■ By equating 
the coefficient of x± ■ ■ ■ x rs in equations (13.90 and (13.71) . we get 

CS = C[ rs ] = ]J Oj + (6 n -6 n _i) Yi °H ^(br s +l-b rs ) J| 0, 

ie[n]\S ie[n-l]\S ie[r s ]\S 

= ]J a i + ( b n -K-l) II «i H h (&r s +l - &r s ), 

ieMV-S 1 ie[n-i]\s 

since aj = TT a» := 1, by definition. Now let S be any nonempty index set. Then 

ie[r s ]\s i€0 



c s = n a * + ( & « ~ n a * h — ( b r s +i - b rs ) n a * 

ie[n]\S ie[n-l]\S ie[rg]\S 

= II a ii II fl i + ( & n-&n-l) II ^( & rs+l-^r s )] 

iGlrslVS 1 i6[n]\[r S ] ie[ra-l]\[r s ] 

= ( II a i) C [r S ] 
ie[r s ]\S 

= C [rs ] J] C N\«- 
ie[rs]\S 

This completes the proof that a nested canalyzing polynomial has to satisfy equation (13.80 . 

Conversely, suppose that q n ] = 1 and equation (13.81) holds for the coefficients of the polynomial 
/ in equation (13.71) . We need to show that / is nested canalyzing. Using Lemma [275| it is enough 
to show that / depends on all n variables and /(of, . . . , %, Xj+i, • • • , x n ) = bj for some 
a,,-, bj G and 1 < j < n. Since c\ n \ = 1, the monomial summand in / and hence 

/ depends on all n variables. Now let 1 < j < n. For any S C [n] such that j S and rs > j, 
we have 

cs = c [rs ] Y[ c [n] \ {i} and c 5u{j} = c [rs] JJ C[n]\{i}- 

ie[r s ]\S ie[r s ]\{SU{j}} 
By pairing with CguU} and with Oru{j} where TC [j — 1], we rewrite the form (13.70 into 

/(a:i, . . . ,x n ) = ( x i c Tu{j} + c T ) Y[xi + (c[ n )\{j} + Xj) c su{j} T[xi (3.H) 

TC|j-l] i€T Sc[n] i&S 

rs>j 

For 1 < j < n, let aj = C[ n ]\{j\. Then 

/(oT, . . . , aJZl, a h x j+1 , ...,x n )= ( c [n]\{j} c Tu{j} + c T ) J| i 1 + C M\«) ( 3 - 12 ) 

Tcy-i] i&T 
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is a constant which we call bj. Hence, by Lemma [2.5) the function / is nested canalyzing. 

Remark 3.4 Observe that the relations in equation A3.8\) leave the coefficients cq and cu, for 
all 1 < % < n, undetermined, as well as the coefficients eg, where S is any of the (n — 1)- 
element subsets of [n] which include n. Furthermore, a Boolean NCF requires that cr n i = 1. 
Since a general Boolean polynomial in n variables has 2 n coefficients, equation 113. 8\) yields 
2 n — 2n equations which have to be satisfied by the coefficients of a Boolean NCF. 

Corollary 3.5 The set of points in F| corresponding to coefficient vectors of nested cana- 
lyzing functions in the variable order xi, . . . ,x n , denoted by V£ , is given by 

^d Cf = U c 0> • • • > °[n]) e F f : C N = 1, c s = c [rs] n c [n] \ {i} , for S C [n]}. (3.13) 

ie[r s ]\S 



The following corollary provides surprisingly simple expressions of the canalyzing input and 
canalyzed output values in terms of the coefficients of the polynomial. 



Corollary 3.6 Let f be a Boolean polynomial given by equation 43.7J) . If the polynomial f is 
a nested canalyzing function in the order xi, X2, • • • , x n , with input values a,j and corresponding 
output values bj,l < j < n, then 



a j= C [n]\{j}, 

b\ = C0 + ciC[ n ]\{i}, 

h 3+i ~ b i = c b'+i] c N\{i+i} + %], 
b n — a n = b n -i + C[ n _x] . 



for 1 < j ' < n 



for 1 < j < n — 1 and 



(3.14) 
(3.15) 
(3.16) 
(3.17) 



PROOF. Equation f)3.14p follows from equation (13.101) . equation (13. 15)) follows directly from 
equation (I3.12p when j = 1. In equation (13.31) . we observe that the variable x n _i appears only 
in the first and second group of products. In particular, q n _i] = — a n + b n — 6 n _i, and hence 
equation (13.171) follows. 

It is left to show (13. 16)) . From equation (13.121) . 



h 3 = i C [n]\{j}CTU{j} + Cr) II ( X + C N\{i}) 

TC\j-l) i€T 

= e ( c N\{i}%i n c n\» + c t) n( i+c N\«) 

TC\j-l) ie[j-l]\T idT 

= C N\{i} C b] £ C [n]\{j} C [j] II C [#)I1( 1 + C [#))+ £ c T Y[(i + c [n] \ {l} ) 

TCy-1] ie[j-l]\T i£T TC[j-l] ieT 

= c H\0> c b1 + J2 °t IK 1 + c N\{i})' 

TC\j-l) ieT 
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since 



e n c n\w ik 1 + c n\w) = n ( c H\{i} + (i + c n\«)) = n 1 

TC\j-i]ie[j-i]\T ieT i€\j-l] ie[j-i] 



Now 



b 3+i - b j = c M\{j+i} c [j+i] + E °t II i 1 + c M\w) ~ c H\{j} c [j] ~ E c t II ( X + C M\«) 

tc[j] ieT rcy-i] ieT 

c N\{i> c b1 + E c t LK 1 + C M\») 

TC[j] ieT 
J6T 



c 



M\{i+i} c [i+i] 



: c N\{i+i} c [i+i] - c M\{j} c [j] + E ( X + c [n]\{j})°[j] II C N\W IK 1 + c M\m) 

Tcy-i] ie[?-l]\T ieT 

: c W\{i+i} c [i+i] ~ c M\{j} c [j] + (1 + c M\{j})c[j] II c M\{i} IK 1 + c N\m) 

Tc[j-i] l e[i-i]\T ieT 

c N\{i+i} c [i+i] _ c N\{i} c [i] + + °[n]\{j}) c [j] 
: c N\{i+i} c [i+i] + %] 



Remark 3.7 Equations l[3.14\ )~ [3-l(>\) imply that the input values ai and output values bi, 
1 < i < n — 1, are determined uniquely by the coefficients of the polynomial f . Also, equation 



( 3.11 ) implies that there are two sets of values for a n and b n which will yield the same nested 
canalyzing function f. Using these facts, we discover Remark \2.9l 

Example 3.8 In Tabled we give some examples of relationships between coefficients of NCFs 
in n variables, nested in the order x\, X2, ■ ■ ■ , x n for some small values of n. 



Table 2 

Some examples of relationships between coefficients of NCFs 



n = 3 


n = 4 


n = 5 


C3 = C13C23C123 
C2 = C23C12 


C4 = C234C134C124C1234 

Cl3 = C134C123 

C24 = C234C124C1234 

C2 = C234C12 


C5 = C2345C1345C1245C1235C12345 
C124 = C1245C1234 
C23 = C2345C123 

C2 = C2345C12 



We now extend Theorem 13.31 to the general case when the variables are nested in any given 
order. For this, we will need to extend the definition of completion of a set S with respect to 
any permutation of its elements. 

Definition 3.9 Let a be a permutation on the elements of the set [n]. We define a new order 
relation < a on the elements of [n] as follows: i < a j if and only if < Let 

S be a nonempty subset of [n], say S = {ii, . . . ,i t } . Let := max{cr~ 1 (ii), . . . , cr _1 (z t )}. 
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The completion of S with respect to the permutation a, denoted by [rg\ a , is the set [rg\ a : = 
{a(l),...,a(r%)}. 

The following corollary is a generalization of Theorem 13 .31 It gives necessary and sufficient 
relations among the coefficients of a NCF whose variables are nested in the order specified by 
a permutation a on [n] . 

Corollary 3.10 Let f G R and let a be a permutation of the set [n]. The polynomial f is a 
nested canalyzing function in the order x a {\), . . . , x a ( n ), with input values a CT (j) and correspond- 
ing output values b a ^), 1 < i < n, if and only if q n ] = 1 and, for any subset S C [n], 

c s = C[rg] CT II C M\M- (3.18) 

w£[r°] a \S 

PROOF. We follow the same argument as in the proof of Theorem 13. 3^ where we impose the 
order relation < a on the elements of [n] and we replace all occurrences of the subscript i by 
a(i) and [r s ] by [rg] CT . 

Corollary 3.11 Let a be a permutation on [n]. The set of points in F^" corresponding to 
nested canalyzing functions in the variable order x a m, . . . ,x a ( n ), denoted by V^ cf , is defined 
by 

K ci = {( c 0> • • • > C N) e Ff : c N = l,c s = c [r g] CT n q n ]\{«}, for S C [n]}. (3.19) 

The following corollary is an extension of Corollary 13.61 and gives the input and output values 
of a Boolean NCF whose variables are nested in the order specified by some permutation a. 

Corollary 3.12 Let f G R and let a be a permutation of the elements of the set [n]. If f is a 
nested canalyzing function in the order x a ^ , . . . , x a r n \ , with input values aj and corresponding 
output values bj, 1 < j < n, then 



a i = c N\Mi)}' for 1 < J < »i - 1 (3.20) 

°l = c + Ca(l)C[n]\{<j(l)}, (3-21) 

~ °j = c b+iU c [n]\W(j+i)} + c [j]a , for 1 < j < n - 1 and (3.22) 

b n - a n = b n -i + C[n-i]„- (3.23) 



PROOF. This follows from Corollary 13.6} where we replace all occurrences of subscript j by 
a(j) and [r] by [r) a . 

Recall that the set V ncf of nested canalyzing functions is the union of the sets V?* of cana- 
lyzing functions with respect to a specified variable order. By Corollaries 13.51 13.111 and the 
correspondence (13. 6p . we have 
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a 



Corollary 13.101 is the starting point for a geometric analysis of the set of all nested canalyzing 
functions. It provides a set of equations that have to be satisfied by the coefficient vectors of 
the polynomial representations of the functions. These coefficient vectors therefore form an 
algebraic variety in the space F^" , which turns out to have very nice properties. In particular, 
it is a so-called toric variety. 



4 Discussion 



Our main contribution in this paper is to connect three different fields of inquiry into Boolean 
functions, which were heretofore apparently unconnected. The equivalence of nested canalyzing 
functions and unate cascade functions relates the electrical engineering point of view of logic 
circuits with the dynamic biological network view, providing a dictionary for results. The 
equivalence of both to a class of polynomial functions brings rich additional mathematical 
structure to the study of both. In particular, the language and concepts of algebraic geometry 
and the rich tool set of computational algebra and algebraic geometry provides a foundation 
that imposes a mathematical structure on the entire class of these functions, which suggests 
an entirely new way of studying them. As an algebraic variety, the class of nested canalyzing 
functions has a very special structure, namely that of a toric variety jij. Toric varieties lie 
at the interface of geometry, algebra, and combinatorics and have a rich structure [2l|]. In 
another paper, we will explore the properties of the toric varieties in the previous section in 
more detail. 



In particular, our motivation for this study originally was the desire to give a characterization 
of nested canalyzing functions as polynomials, which could be used as part of the model 



selection algorithm in [111 ]. That is, we are interested in giving an efficient criterion which 
allows our symbolic computation algorithm to preferentially pick nested canalyzing functions 
rather than general polynomials. The characterization of this class as a toric variety is the 
first important step in this direction. 



It deserves mention that the connection to unate cascade functions was discovered in a round- 
about way. We first established the parametrization of nested canalyzing functions by special 
polynomials. The structure of these polynomials makes it easy to count how many there 
are for a given number of variables. After carrying out this counting procedure for the first 
few numbers resulted in a sequence of integers which we submitted to N. Sloane's integer se- 
quence database (http://www.research.att.com/~njas/sequences/). One of the matching 
sequences was that for the number of unate cascade functions. 
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